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EXPERIMENTS  ON  LOW  ASPECT  RATIO  HYDROPLANES 
TO  MEASURE  LIFT  UNDER  STATIC  AND  DYNAMIC  CnNDTTTnNS 

By  B  Ward  and  A  R  J  M  Lloyd 

1.  INTRODUCTION 

In  support  of  SUBSIM  (Reference  1)  the  computer  program  which  simulates 
submarine  manoeuvres,  there  was  a  need  for  data  on  dynamic  effects  on  lift 
forces  for  low  aspect  ratio  hydroplanes.  A  thorough  literature  search  was 
carried  out;  previous  papers  usually  described  two  dimensional  experiments 
at  much  higher  frequencies  than  those  of  relevant  interest.  It  was 
therefore  decided  to  perform  some  dedicated  experiments  to  obtain  the 
required  data. 

2.  EXPERIMENT 

Three  NACA  0020  hydroplanes  were  used  of  aspect  ratios  1,  1.5  and  2 
(Figure  1).  Each  hydroplane  had  a  chord  length  of  0.26  m  and  span  of 
0.13  m,  0.195  m  and  0.26  m  respectively.  A  strain  gauged  stock  at  the 
quarter  chord  position  of  the  hydroplane  was  attached  to  a  servo  to 
oscillate  the  hydroplane  (Figure  2).  The  same  stock  was  used  for  all  three 
hydroplanes.  Strain  gauge  bridges  on  the  stock  were  calibrated  to  give 
normal  and  tangential  forces. 

The  experiment  was  conducted  in  the  Circulating  Water  Channel  (Figures  3 
and  4).  The  servo  was  attached  below  the  CWC  floor  and  the  position  of  the 
floor  was  altered  to  accommodate  each  hydroplane.  A  1  mm  space  was  left 
between  the  floor  and  the  hydroplane  root.  Horizontal  flow  velocity  was 
maintained  at  2.5  m/s. 

3-  DATA  ACQUISITION  AND  DATA  ANALYSIS 

The  DATS  computer  package  was  used  for  both  data  acquisition  and  data 
analysis.  DATS  consists  of  a  number  of  modules  which  enable  mathematical 
functions  to  be  applied  to  files  of  time  history  data  (see  DATS  User 
Manual,  Reference  2). 

Velocity  of  the  flow  was  measured  upstream  of  the  hydroplane  using  laser 
doppler  velociraetry  techniques. 

Velocities  from  the  laser,  achieved  angle  from  the  servo,  and  normal  and 
tangential  forces  from  the  stock  were  acquired  simultaneously.  The  lift 
was  derived  from  (see  Figure  5) 

L  =  N  cos  a  +  T  sin  a  (1) 

The  upper  tangential  force  gauge  was  damaged  when  the  experiment  was  being 
rigged  and  no  measurements  of  total  tangential  force  were  possible.  The 
readings  from  the  lower  gauge  were  doubled  to  allow  equation  (1)  to  be  used 
to  calculate  the  lift  force.  This  procedure  was  not  strictly  correct  but 
will  not  have  introduced  any  large  errors  because  the  contribution  of  the 
tangential  force  to  the  lift  is  very  small. 
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(2) 


Lift  coefficients  were  calculated  using  the  instantaneous  velocity. 

C  -  — ^  - 
L  "  -ipU^A 

4.  STATIC  RUNS 

Mean  values  of  a  and  were  taken  for  static  runs  over  the  acquisition 
time;  typically  50  seconds  (Figure  6), 

Figures  7  to  9  show  plots  of  non  dimensional  lift  against  angle  of 
attack  for  aspect  ratios  1,  1.5  and  2  respectively. 

Figure  10  shows  a  plot  of  stall  angle  against  aspect  ratio. 

5.  OSCILLATING  RUNS 

Oscillating  runs  (Figure  11)  were  carried  out  at  frequencies  from  0.055  Hz 

to  0.386  Hz  corresponding  to  ^  values  from  0.036  to  0.252. 

Fre^encies  were  chosen  to  correspond  to  full  scale  angular  rates  in  the 
region  of  5  degs/sec.  The  frequency  was  limited  to  ensure  that  rate 
limiting  did  not  occur  in  the  servo.  The  average  C.  amplitudes  were 

ij 

calculated.  Plots  are  shown  in  Figures  12  to  14.  Little  variation  with 
frequency  was  apparent. 


Figures  15  to  17  show  plots  of  phase  angles  against  ^  for  three  different 

aspect  ratios.  The  phase  angles  are  positive  indicating  that  lift  leads 
incidence.  Note  the  linear  trend  in  the  data. 

Figure  18(i)  to  18(iii)  shows  continuous  plots  of  C^^  against  incidence  for 

three  oscillating  runs  with  aspect  ratio  2  at  different  frequencies. 
Incidence  amplitude  is  35.5  degrees  in  all  cases  (above  the  static  stall 
angle).  Note  at  the  lower  frequency  stall  is  evident.  At  the  higher 
frequencies  the  lift  coefficient  is  higher  at  zero  angle  corresponding  to 
greater  phase  lead.  Stalling  does  not  seem  to  occur  at  the  high 
frequencies. 

6.  RAMPING  MOTIONS 

Figure  19  shows  the  results  of  three  ramping  motion  experiments  to  an 

incidence  of  41  degrees  at  6  degs/sec  (^  =  0.011)  for  an  aspect  ratio  of 

1.5.  In  each  case  the  lift  coefficient  follows  the  angle  of  incidence  to 
the  maximum  static  value.  In  other  words  the  hydroplane  does 
not  immediately  stall.  After  some  time  flow  separation  occurs,  the 
hydroplane  stalls  and  the  lift  collapses  to  the  value  given  in  Figure  8. 
High  lift  is  maintained  for  an  apparently  arbitrary  time. 
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7 • 1 •  Representation  of  by  Second-Order  Linear  Equations 


It  is  assumed  that  under  dynamic  conditions  the  equation  for  lift 
coefficient  is 


«  ^  i_  •Cl  ••  c 


(3) 


where  a^a  is  the  static  lift  curve  slope,  a^a  is  the  effect  of  angular  rate 


and  ^20  is  the  acceleration  effect. 


The  following  equations  can  then  be  derived 


/  c_  2  c  .  V 

“  0  ^  2  2  “  ^  ~  ^  1  u“  ^  ^  0  ^  L 


where  is  the  a  amplitude 
and  is  the  amplitude 

l-Q  L 

and 


«  =  tan 


-1 


liJL 


2  c‘ 

■  ®2"  JJ2 


(4) 


(5) 


Using  the  experiment  results  and  the  least  squares  plane  method  the 
coefficients  in  equation  (4)  were  derived  for  each  of  the  three  aspect 
ratios; 


a 

e 

^2 

^0 

1 

0.37 

1.12 

1.72 

1 .5 

0.44 

1.3 

1.95 

2 

0.47 

1.4 

2.12 

The  experiment  gave  lower  values  than  the  empirical  equation  (6)  derived  by 
Whicker  and  Fehlner  (Reference  3)  at  aspect  ratio  1.5  and  2: 


(6) 


1 .8  IT  a 


^0  = 


1.8  + 


per  radian 


That  particular  experiment  used  a  NACA  0015  section  as  opposed  to  the 
NACA  0020  used  in  this  study.  NACA  0020  sections  give  reduced  values  of 

compared  with  NACA  0015  (Reference  4).  However  at  aspect  ratio  1  the 
Whicker  and  Fehlner  equation  gives  lower  values  than  this  experiment 
justifying  the  conclusion  made  in  Reference  5  -  'that  the  lift  curve  slope 
at  very  low  aspect  ratios  is  seriously  under  estimated  by  equation  (6)  . 
Although  that  conclusion  was  based  on  only  two  results,  the  results  from 
this  experiment  give  added  weight  to  the  statement.  Using  these 
assumptions  and  taking  the  lift  curve  slope  to  be  zero  at  zero  aspect  ratio 
formed  the  basis  under  which  a  more  general  equation  was  derived  for 
coefficient  ag  ie  the  lift  curve  slope: 

ag=  163  tanh  (0.5  a^+  2.4)  -  160  per  radian  a^s  2 

(7) 

_ ^  ^e _  per  radian  a^  >  2 

/  2  4 

1.8  +  cos  074+  a^/cos  O 


A  plot  can  be  seen  in  Figure  20, 

As  can  be  seen  from  Figures  7  to  9  the  equation  fits  the  linear  part  of  the 
experiment  data  for  static  conditions  very  well. 

Equations  were  derived  to  fit  the  coefficients  a^  and  a2  in  equation  (3). 

Only  three  points  are  available  and  it  was  assumed  that  there  are  no 
dynamic  effects  if  aspect  ratio  is  zero.  It  was  assumed  that  the  function 
approaches  limits  determined  by  the  trend  of  the  data  available.  A  plot  of 
the  equations  for  coefficients  and  a2  can  be  seen  in  Figures  21  and  22 
respectively. 

a^  =  99.3  tanh  (0.61  a^  +  2.42)  -  97.7  (8) 

a2  =  1.21  tanh  (0.78  a^  +  0.68)  -  0.71  (9) 

The  satisfactory  fit  of  the  equations  can  be  seen  in  Figures  12  to  17. 

Equation  (5)  for  phase,  compared  with  experimental  results  in  Figures  15  to 
17,  is  best  for  a^  =  2  and  adequate  for  the  lower  aspect  ratios. 


7.2.  Equation  for  Stall  Angle 

An  empirical  equation  was  derived  from  the  experiment  data  on  stall  angle 


(10) 


2 

o_„  =0.113  a  -0.615  a  +1.32  radians 
ST  e  e 

This  gives  higher  values  for  all  three  aspect  ratios  than  that  calculated 
by  the  formula  used  in  SUBSIM  (Reference  1) 

a_„=  0.075  a^  -  0.445  a  +1.05  radians  (11) 

ST  e  e 


which  is  derived  from  the  data  in  Reference  3.  These  equations  are 
compared  with  the  experiment  data  in  Figure  10. 

The  above  equation  is  valid  for  1  <  a^  <  2. 

8.  EFFECTS  DUE  TO  THICKNESS  CHORD  RATIO 
8.1.  Correction  for  Lift  Curve  Slope 

From  results  quoted  in  Reference  4,  (Figure  23)  it  can  be  seen  that  a 
change  in  thickness  ratio  has  an  effect  on  lift  curve  slope.  From  these 
data  the  following  formula  has  been  derived  to  correct  for  thickness/chord 
ratio. 


.Q=  (0.96  +  1.7  t/c  -  7.1 


^2/  2, 
t  /  c  ) 


(12) 


The  formula  is  valid  for  data  in  the  range  t/c  =  0.09  to  0.35. 

8.2.  Correction  for  Stalling  Angle 

Reference  3  gives  stall  angle  data  for  square  tip  fins.  Using  those 
results  and  the  results  from  this  experiment  a  general  equation  for  stall 
angle  was  derived  correcting  for  thickness  ratio. 

=  (0.36  +  3.2  t/c)  (13) 


9.  EFFECTS  DUE  TO  FAIRED  TIPS 

An  unpublished  report  by  the  Stevens  Institute  of  Technology  discusses  the 
fact  that  square  tip  hydroplanes  give  greater  lift  forces  than  faired  tip 
hydroplanes.  This  was  also  noted  in  Reference  3  by  Whicker  and  Fehlner. 
The  conclusion  reached  in  the  Stevens  report  was  that  the  hydroplane  tip 
outboard  of  the  trailing  edge  span  can  be  assumed  to  carry  no  lift.  This 
should  be  taken  into  account  in  SUBSIM  by  specifying  the  hydroplane 
outreach  to  exclude  any  faired  tip. 
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10.  RANDOM  SIGNAL  EXPERIMENTS 


Following  the  experiments  described  in  this  report  a  further  experiment 
(Reference  6)  was  conducted.  The  1.5  aspect  ratio  hydroplane  was 
oscillated  using  a  random  signal.  Figure  24  shows  a  plot  comparing  the 
lift  coefficient  derived  from  equation  (3)  and  a  time  history  of  from 

the  random  signal  experiment.  The  comparison  is  very  good,  confirming  the 
validity  of  equation  (3). 

11 .  CONCLUSIONS 

The  experiment  to  examine  dynamic  effects  on  oscillating  hydroplanes  showed 
the  effects  to  be  small  but  worth  noting.  The  most  interesting  aspect  of 
the  experiment  was  the  delayed  stall  noted  during  the  ramping  motion 
experiments . 

The  rapid  change  in  flow  velocity  at  the  surface  of  the  hydroplane  as  it  is 
ramped  over  causes  a  rapid  change  in  pressure  gradient.  There  is  a  delay 
in  the  build  up  of  the  adverse  pressure  gradient  which  allows  the 
hydroplane  to  support  greater  lift  than  during  the  steady  state.  This  was 
noted  by  Ericsson  in  Reference  7.  According  to  Moore  (Reference  8)  this 
delay  allows  the  hydroplane  to  respond  to  the  change  of  angle  of  attack 
without  stall  which  contributes  to  the  lift  overshoot. 

Above  the  stall  angle  oscillating  motions  show  increased  values  for  Cj^ 
compared  with  static  results  due  to  dynamic  effects. 

No  significant  variation  in  was  noted  over  the  range  of  frequencies 
examined.  Phase  angles  are  small;  a  phase  lead  of  approximately  10  degrees 
occurs  at  the  highest  frequency  of  0.386Hz  (^  =  0.252).  This  may  be 
significant  for  autopilot  design. 

Effects  due  to  thickness  ratio  and  faired  tips  should  be  considered  when 
incorporating  any  empirical  formula  derived  from  this  experiment  into 
SUBSIM. 

Empirical  formulae  fit  experimental  results  very  well.  These  formulae 
should  be  incorporated  into  the  SUBSIM  mathematical  model  and  assessed  by 
further  validation  of  submarine  predictions  against  full  scale  trial 
results . 
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FIG.1  THREE  NACA  0020  HYDROPLANES 


FIG. 3  EXPERIMENT  IN  CIRCULATING  WATER  CHANNEL 
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FIG. 4  EXPERIMENT  IN  CIRCULATING  WATER  CHANNEL 
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FIG.  9.  LIFT  COEFFICIENT- ASPECT  RATIO  2,  STATIC 

CONDIT  ION 


FIG-  lO.  STALL  ANGLE  -STATIC  CONDITION 
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FIG.  12.  LIFT  COEFFICIENT  ASPECT  RATIO  I,  DYNAMIC 
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FIG. 13,  LIFT  COEFFICIENT- ASPECT  RATIO  1-5,  DYNAMIC 

CON04TION 
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FIG.  15.  PHASE  ANGLE  -ASPECT  RATIO  1 
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FIG.  16.  PHASE  ANGLE -ASPECT  RATIO  1-5 
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FIG.  17.  PHASE  ANGLE-  ASPECT  RATIO  2 
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FIG.18(i)  DYNAMIC  CONDITION  -  ASPECT  RATIO  2 
35.5  DEGREES  ,  FREQUENCY  0.055HZ 
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FIG.  20.  EMPIRICAL  FORMULA  FOR  COEFFICIENT  Qq 
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